Macroscopically distinct quantum superposition states as a bosonic code for 

amplitude damping 
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We show how macroscopically distinct quantum superposition states (Schrodinger cat states) may 
be used as logical qubit encodings for the correction of spontaneous emission errors. Spontaneous 
emission causes a bit flip error which is easily corrected by a standard error correction circuit. The 
method works arbitrarily well as the distance between the amplitudes of the superposed coherent 
states increases. 
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I. INTRODUCTION 



In quantum information theory logical states are en- 
coded as two orthogonal pure states 0. The simplest 
example is provided by a single two level system. The 
ground state \g) and excited state |e) can then encode 
logical and logical 1 respectively. The ability to form 
a coherent superposition of logical states is why we refer 
to logical states as qubits rather than simply as bits Q . 

Quantum computation gains its power through the po- 
tential ability to unitarily manipulate a coherent super- 
position of large collections of physical systems each en- 
coding a sing le qubit @. There is no fundamental reason 
to restrict oneself to physical systems with two dimen- 
sional Hilbert spaces for the encoding. It may be more 
natural in some contexts to encode logical states as a su- 
perposition over a large number of basis states. When 
the system supporting the qubit encoding is coupled to a 
perturbing environment an extra unwanted, and possibly 
uncontrollable, unitary interaction is introduced which 
can appear as an error in the encoded information. The 
coupling to the logical basis determines the type of logi- 
cal error. While the coupling to the environment is fixed, 
we are free to choose how we encode the qubits, hence 
the choice of basis for the logical encoding may change 
the kind of error introduced. For example, with a single 
qubit, a bit flip in one logical encoding basis can appear 
as a phase flip in another B. This is relevant, as some 
kinds of errors are easier to fix than others. Chuang and 
Yamamoto || recently introduced a qubit coding for two 
bosonic modes. 

These modes could be two optical modes or two vi- 
brational modes of a single trapped ion. A particularly 
difficult source of error for bosonic modes arises from 
exponential decay of the energy. In a single mode, for 
example, one could use the ground state and first excited 
state as the logical basis. While the ground state is in- 
variant under decay, the first excited state will 'reset' to 
the ground state in a single decay event. Such an error 
can in general be corrected by a five qubit code || . How- 
ever the code of Chuang and Yamamoto enables a more 
efficient error correction. 



In this paper we give an example of how a careful choice 
of the coding scheme can make a difficult error correction 
task simpler. Our example is based on a quantum code 
for a single bosonic mode that enables amplitude damp- 
ing or amplification to be corrected as a bit flip error. 
The code is based on quantum superpositions of bosonic 
coherent states, the so called 'cat states' 0. Our coding 
scheme is not exact for very small amplitude coherent 
states, but improves exponentially when amplitudes are 
greater than unity. We demonstrate a completely uni- 
tary, adiabatic method to generate the cat states of our 
coding scheme. 



II. CAT STATE ENCODING FOR AMPLITUDE 
DAMPING 

Let | a) be a coherent state for a single bosonic de- 
gree of freedom. We then define two orthogonal states as 
symmetric and anti-symmetric superposition of coherent 
states by 



\S)=Af+(\a) + \-a)) 
\A)=Af-(\a)-\-a)) 



(la) 
(lb) 



where a is an arbitrary complex number. The normali- 
sation constants are given by 



A4 



-1/2 



(2) 



It is easy to verify that the symmetric cat state, \S) con- 
tains only the even energy eigenstates, while the anti- 
symmetric cat state \A) contains only the odd energy 
eigenstates. This feature is independent of a. The two 
states are orthogonal and we are led to the following log- 
ical encoding for a single qubit, 



|0)l = |S) 
\1)l = \A) 



(3a) 
(3b) 



Under free dynamics, the coherent state evolves as 
\a(0)e~ ZUJt ) , however the two cat states remain orthogo- 
nal and thus the logical encoding of the qubit is invariant 
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under free dynamics. Therefore we can transform to the 
interaction picture rotating at frequency oj. 

The amplitude damping model is the standard one for 
a bosonic mode, of frequency u>, weakly coupled to a 
zero temperature heat bath j5|^] in the Born and Markov 
approximation. The system obeys the following master 
equation in the interaction picture 



d P 7/ t t t > 

— = — (Zapa 1 — a 1 ap — pa' a) 

The solution to this equation may be written as 



P (t) = J2r k (t)p(o)rt(t) 



(4) 



(5) 



with 



T*(t) = £ J (t) (V(t)) (n - k)/2 (l - V(t)) k/2 \n k)(n\ 

n— k V ^ ^ 

(6) 

and rj(t) = e~ 7 * is the probability that the state is unde- 
cayed up to time t. 

Our objective is to correct for at most one decay event 
over some characteristic time. In which case we only 
need to consider the two terms corresponding to To and 
Ti. Coherent states remain coherent under amplitude 
damping, and in particular we have that, 



(1-'))M/2| 

Ti|a) =a^e-( 1 -">HV2|^ a) . 



(7a) 
(7b) 



It is then easy to see that a single decay event will cause 
an even cat state to flip to an odd cat state and vice 
versa. It is this feature that we are attempting to exploit 
through our code states, so that a single decay event 
will correspond to a bit flip. A no decay event essentially 
leaves the state unchanged. These statements are strictly 
only true for cats with an infinitely large coherent ampli- 
tude, however we now show that only small amplitudes 
are sufficient for practical purposes. 

An error correction code must satisfy the following con- 
ditions fol 



L (p\TlT l \q) L = fovp^qork^l 

L (p\r{r k \ P ) L = p k for P = o,i 



(8a) 
(8b) 



where p, q are or 1 , and P k is a constant that depends 
only on k. The first equation requires that all erroneous 
states are orthogonal and the second requires the proba- 
bility for each event (no decay or one decay) to occur to 
be independent of the logical state. It is easy to see that 
the cat state encoding satisfies the first condition. The 
second condition however requires more careful consider- 
ation. Using the conditional states given above we find 
that 



3 — 2r]a 2 



l (0|T Tq|0) l = 1 + e- 2 "" 2 
L (1|T T |1) L 1 
L<0|Tfri|0) L 
L (l|T t 1 T 1 |l) L 



1 



-2r/Q 2 



(9a) 
(9b) 



Each of these ratios should ideally be unity, but the de- 
parture from ideality is insignificant even for such a small 
value as a — 3. For example, with r\ = 0.9 we find ( |9a| ) 
gives 1.00149 for a = 2 but for a = 3 it gives 1.000000184. 
While @ gives 0.9985079 for a = 2 and 0.999999815 for 
a = 3. If we increase the amplitude to a = 5 the depar- 
ture from ideality is undetectable. Therefore the logical 
qubits are encoded in a manner that enables amplitude 
decay to be corrected to any desired degree of precision. 

We can see that after many spontaneous emission 
events the amplitude will eventually decay away to zero. 
If the coherent amplitude is too small then the ratios 
( pa| ) and (9b) will deviate significantly from unity. It is 
therefore necessary to have a sufficient initial amplitude 
to allow computation for a reasonable amount of time 
and to know when it is prudent to reset the states. 

It is possible to determine the time scale over which 
the states will be useful by considering the ratio (|9^). 
This ratio should not be significantly different from 1 for 
the encoding to work, so we allow the difference to be no 
greater than a small tolerance. The term responsible for 
any deviation of the ratio is exp(— 2-ija 2 ) which we desire 
to be small enough such that the ratio is within toler- 
ance. This implies that r]a 2 has to be greater than some 
limiting value determined by the tolerance, below which 
the state must be reset. Therefore, given a certain error 
rate, initial coherent amplitude and desired tolerance we 
have sufficient information to calculate the time available 
for computation before reset. 



III. LOGICAL OPERATIONS ON CAT STATES 

A logical encoding is useless if we cannot implement 
one and two qubit operations on the encoded states. We 
now show how this can be done for the cat state encoding 
defined above. The particular form of qubit operations 
depends upon the particular physical realisation of the 
bosonic mode. For the purposes of illustration we sim- 
ply postulate particular bosonic interactions to achieve 
the required gate operations. We will show that the 
Hadamard transform may be implemented by simple dis- 
placement of a single bosonic mode, while the two qubit 
operation may be realised by a mutual phase shift inter- 
action term which commutes with the number operator 
of each bosonic mode. 

If the bosonic mode is subject to a classical driving 
force the Hamiltonian describing this process in the in- 
teraction picture is 



H D = h(p<f + (3* a) 



(10) 
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where (3 is the complex amplitude of the driving force. 
Let us now choose (3 as real (in general we choose (3 to 
be 7r/2 out of phase with a). 

For a given cat state amplitude we can choose the driv- 
ing amplitude such that 



9 = af3t 



(11) 



where t is the length of time the driving force is applied. 

If the even cat state (encoding \Q)l) is driven we find 
that 



-iH D t/h 



\0) L = cos0|O) L -zsin0|l) L 



(121 



A displacement of this kind shifts the "cat" very 
slightly by an amount (3 in a direction orthogonal to the 
orientation of the cat state in phase space. The transfor- 
mation is approximately equivalent to a Hadamard trans- 
form of the single logical qubit when 8 = 7r/4 (in the limit 
of large a and small (3) and will suffice as a universal one- 
qubit gate. We will refer to this as a H-gate. 

The simplest way to realise a two qubit universal gate 
is via the two mode interaction Hamiltonian 



Hp = hxa'ab^b 



(13) 



where a, b represent the mode amplitude operators for 
the two bosonic modes of interest. We choose the inter- 
action time t such that %i = ir. As the |0)l only has 
even bosonic number while only has odd bosonic 
number, we find that the interaction leaves the states, 
|0)La|0)i b , |0) La |l) L fc, |l)z, a |0)i6 unchanged, but the 
state in which both modes encode a |1)l transforms as 



„— iita* ab^b | -i \ | 
e 1/ia 



1>« = -|1> 



La|l)Lfc 



(14) 



This kind of conditional phase shift operation suffices for 
a universal two qubit gate. We will refer to this as a 
P-gate. 

Using the one and two qubit gates described above, we 
can construct a controlled-not (CN) gate. Let mode a 
code the control bit and mode b code the target bit. A 
CN gate is then made by applying a H-gate to the target, 
then coupling the target and the control by a P-gate and 
finally applying another H-gate to the target. 

We have shown that simple one mode and two mode 
transformations may be used to construct universal com- 
putational gates for a cat state logical encoding of bosonic 
systems. Amplitude damping appears as a simple bit-flip 
error in this encoding and thus a three qubit code can be 
used to correct it. This leads to relatively simple fault tol- 
erant implementations of the gate operations described 
above using three coupled bosonic modes. 



IV. UNITARY CONSTRUCTION OF CAT STATE 
ENCODING 

The cat state encoding described in this paper will be 
of little use if we cannot encode our logical bits by uni- 
tary transformations. Unfortunately all previous schemes 



to generate cat states are based on an entanglement be- 
tween a bosonic mode and a two level atom and require 
a measurement readout The cat state produced is 

conditional on the two, mutually exclusive, results of this 
measurement, and we are equally likely to get an even cat 
as an odd cat. This method of encoding would randomly 
assign logical bits and is of little practical use. We now 
describe a unitary, although adiabatic, method to gener- 
ate the two kinds of cat state used to encode the qubits. 
Consider the Hamiltonian. 



H 



NL 



hx(a 



(15) 



which could describe a Kerr nonlinearity for an optical 
bosonic mode or the self-interaction of a single trapped 
ion driven at the carrier frequency [JL3| , in which case % 
is proportional to the fourth power of the Lamb-Dicke 
parameter. The Hamiltonian in (15) has two degener- 
ate ground states which are the ground state, |0) and 
first excited state, |1) of a single bosonic mode. In both 
cases the eigenvalue is zero, however each of these ground 
states is distinguished by the parity operator, where |0) is 
even and |1) is odd. If we now consider the Hamiltonian 



He 



H 



NL 



h K (a 2 + (a^) 2 ) 



(16) 



with k > 0. Noting that the cat states |0)l, are 
eigenstates of a 2 , it is easy to see that these sa me st ates 
are degenerate eigenstates of He when a = \[k]x and 
the eigenvalue is —Tlk 2 /x- While the cat states are degen- 
erate eigenstates of He they are distinguished by their 
parity. The adiabatic theorem now enables us to predict 
that the even or odd initial eigenstates of -Htvl, |0), |1), 
will evolve respectively into the even or odd eigenstates, 
|0)ij |1)l of He as we slowly turn on k from zero to a 
final target value. Thus we have a unitary method to 
code either a logical zero or logical one as a cat state 
by choosing to start from a bosonic ground state or a 
bosonic first excited state. 

The adiabatic theorem is exact only in the case of infi- 
nite slowness, which is of little use for logical encoding in 
quantum computation, so what matters is how well we 
can do in practise. To test this we consider two different 
ways to vary n in time; linear and nonlinear. 

The linear variation considered here consists of sim- 
ply increasing k according to k = t. The function 
k = fco tanh 2 (Xt) was used in the nonlinear case due to 
the advantageous shape of the tanh 2 function. 

Figure |] illustrates the fidelity versus time variation 
of the state starting from the |1) Fock state with re- 
spect to the equivalent cat state of mean photon number 
a = \J K /x- The fidelity is measured as the modulus 
squared of the dot product of the evolving state with the 
cat state. The notable features of Figure [j] are the fi- 
delity oscillations, the "steady-state" fidelity and the rel- 
ative characteristics of the linear and nonlinear methods 
of varying n. 
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The fidelity oscillations result from carrying out the 
adiabatic evolution faster than as required for exactness 
by the adiabatic theorem. As the system evolves from 
t = the fidelity will tend away from unity as the state 
\ip) evolves away from the relevant cat state. Continued 
evolution eventually causes to more closely resem- 
ble the equivalent cat state with the fidelity increasing 
accordingly. The retreat and approach of the evolving 
state with respect to the cat state causes the oscillations 
seen in Figure |[ 

The oscillations are damped by k until a "steady state" 
is reached with constant fidelity. The "steady state" fi- 
delity is determined by how quickly k is increased from 
t = 0; a slower initial increase implies a greater final fi- 
delity. Hence there arc two effects occurring with Kj\ as k 
increases, oscillations in fidelity are suppressed, and the 
faster k is increased initially, the lower the steady state 
fidelity of \ip). 

For this adiabatic process to be useful we have two 
aims; a steady state in a reasonable amount of time, and 
a state \if>) as close as possible to the desired final state. 
For linear k these are complementary, although for non- 
linear k we can choose a function which can achieve both 
aims, hence the use of the tanh 2 function (s-curve). The 
s-curve has the properties that it starts slowly, thus giv- 
ing a high final fidelity, and later damps the system very 
quickly to give a useful final state in a reasonable amount 
of time. If the variables fco an d A are chosen carefully 
then it is possible to obtain a fidelity of almost unity in 
a usefully short time. 

We thus conclude that the unitary logical encoding in 
terms of cat states may be performed with almost arbi- 
trary accuracy using this adiabatic method. 



at an antinode of an optical standing wave tuned to the 
atomic frequency; the carrier frequency. In an interac- 
tion picture at frequency v the interaction Hamiltonian 
is, 

fin 4 

H I = -hfln 2 a^aa x + h—^(a^) 2 a 2 a x (17) 

where Q is the Rabi frequency and r\ is the Lamb-Dickc 
parameter. The first term in this expression suffices to 
build a CN gate between the cat state and the electronic 
state. If we choose the interaction time appropriately we 
can apply the transformation 

U — exp(—iwa' aa x ) (18) 

When this acts on an even cat state it corresponds to the 
identity on the electronic system. When it acts on an odd 
cat state it corresponds to a 7r pulse in the electronic sys- 
tem. If we code our electronic qubits as \g)i —* |0)j and 
|e)i — > The unitary interaction in Jig ) will effect 

a CN gate with the bosonic mode acting as the control 
and the electronic mode acting as the target. Thus joint 
excitation on the carrier frequency of the two ion sys- 
tem will produce the double CN gate in the first part of 
Figure ||. The final double CN gate in which the vibra- 
tional mode becomes the target can easily be produced 
with the same Hamiltonian with H-gates either side. This 
procedure would enable a cat state, once produced, to be 
protected from single decay events. 



ACKNOWLEDGMENTS 



V. DISCUSSION AND CONCLUSION 

We have shown that the even and odd cats states may 
be used as a robust qubit encoding for a single bosonic 
mode subject to amplitude damping. A single decay 
event will then appear as a simple bit flip error. We 
have also shown how the states may be prepared unitar- 
ily and how one qubit and two qubit universal quantum 
gates may be realised. We now turn to an assessment of 
how practical the scheme is for present technology. To 
be specific we will consider the case in which the bosonic 
mode is the centre-of-mass vibrational state of a single 
trapped ion. Cat states have been produced in these sys- 
tems using a conditional measurement scheme |llj . 

Given a cat state it is straightforward to protect it 
against decay using two additional qubits. These could 
be the electronic states of two ions in the trap. The 
error correction circuit for a bit flip is well known and 
is given in Figure To implement the gate we need 
to implement a CN gate between the vibrational state 
and the electronic states of the two ions. Following de 
Matos Filho and Vogel |13 we consider an ion trapped 
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FIG. 1. Fidelity-time evolution, tj) starting from |1). 
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FIG. 2. 3-qubit circuit to correct bit-flip errors. 
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